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A tractable A''-state Rabi Hamihonian is introduced by extending the parity symmetry of the 
two-state model. The single-mode case provides a few-parameter description of a novel class of 
periodic systems, predicting that the ground state of certain four-state atom-cavity systems will 
undergo parity change at strong coupling. A group-theoretical treatment provides physical insight 
into dynamics and a modified rotating wave approximation obtains accurate analytical energies. 
The dissipative case can be applied to study excitation energy transfer in molecular rings or chains. 
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Interactions between spin systems and harmonic oscil- 
lators (boson modes) have been studied for over 70 years 
[ll-Q. One of the most well-known, the quantum Rabi 
model [l[, is a phenomenological Hamiltonian describ- 
ing interactions between a two-level system and a cavity 
mode. The model has also formed a basis of understand- 
ing for exciton-phonon interactions Q and, along with 
its multi-mode extension, has numerous established ap- 
plications in chemistry and physics (see 5|-|7[ and refs. 
therein). The Jaynes-Cummings (J-C) |3| model is ob- 
tained by taking the Rabi model in the rotating wave ap- 
proximation (RWA), where the "counter- rotating" terms 
are ignored (see e.g. 3). While the J-C model is suffi- 
cient to study small atom- field coupling, the RWA breaks 
down at large coupling [13] and the full Rabi model 
is needed. Experimental techniques have accessed these 
strong-coupling regimes and there is much ongoing 
interest in future experimental realizations in both cavity 
111 and circuit ^ QED. 



Many-site spin-boson interaction, e.g. multi-state 
atom-cavity interaction |14| or excitation energy trans- 
fer in multi-chromophoric systems [isj, continues to be 
a subject of significant interest, dictating a need for 
extensions of the two-state model. Extensions of the 
J-C model have been studied extensively 16-1^, but 
are no longer applicable in the strong-coupling regime. 
Exciton-phonon generalizations which extend the par- 
ity/reflection symmetry of the Rabi system j^^] are nei- 
ther tractable nor applicable to atom-cavity systems. 
Most importantly, the Rabi model is the single-mode 
version of a dissipative (infinite-mode) spin-boson model 



21j . signifying that light- matter interaction is a sim- 
plified manifestation of a more fundamental interaction 
between a two-state system and a dissipative environ- 



ment. Previous dissipative |22h25| generalizations have 



neither extended the symmetry nor preserved this cor- 
respondence. Motivated by these properties, this Let- 
ter presents a symmetry-preserving A'^-state extension of 
the Rabi model. The extension includes counter-rotating 
terms in a rigorous, intuitive, and mathematically man- 



ageable way, using a minimal number of parameters and 
paving the way for applications to multi-level atom-cavity 
experiments at both weak and strong coupling. A group- 
theoretical approach 0] provides numerical advantages 
and physical insight into dynamics of the single-mode 
case. The symmetric generalized RWA 0] is applied to 
obtain accurate analytical energies/eigenstates valid for 
strong coupling. The above procedures are significantly 
simplified via the generalized spin matrices |26|, provid- 
ing a new tool for the treatment of general periodic sys- 
tems. The corresponding infinite-mode extension can in 
turn be applied to periodic dissipative A^-state systems. 




FIG. 1. (color online) (a) The partially diagonalized (left 
panel), position (center), and energy (right) representations 
form three equivalent interpretations of the two-state Rabi 
model. (b),(c) analogously depict the Ai'-state Rabi model 
from Eq. ((6]) with N = 3, 4, respectively. The parame- 
ters J, K form the energies of the system, {p, 5, c} are each 
system's conserved quantum numbers, and //jv,n are boson 
chains (defined in text). 

This work will present the Af-state Rabi model's physi- 
cal motivation in two different representations, discussing 
the A^ = 3,4 cases from the viewpoint of atom-cavity 
physics. A discussion of the conserved quantum numbers, 
a third symmetry-based representation, and dynamical 
properties of the single-mode case follows. The remain- 
ing space is devoted to a brief description/application of 
the dissipative case. For reference, the three representa- 
tions discussed here are graphically outlined in Fig. [l]for 
A^ = 2,3,4. 



Position representation. — Consider the Hamiltonian 
H = Hgys + Hfieid + Hint describing an TV-state sys- 
tem interacting with electric and magnetic fields (E and 
B) of a cavity mode of frequency to and wavenumber 
The Hamiltonian will be formulated in the position rep- 
resentation, where the interaction Hi„t is on the (spin) 
diagonal. This differs from the traditional introduction 
of the Rabi model in the energy representation, where 
the interaction is off-diagonal (pg. 194 of [I])- While the 
two formulations are equivalent for the two-state case, 
this version provides a natural symmetry-based exten- 
sion. Assume a class of systems defined by 

JV-l 

Hsys = E E -^^^ + ^1 + \^ + k){n\) , (1) 

k n=0 

where {\n)}^^Q form a complete set of position eigen- 
states and k sums over all neighboring sites. Transform- 
ing H into the energy representation, i.e., diagonalizing 
Hsys , would obtain an A''-state "atom" with energies de- 
termined by the parameters J^. The interaction is 

H^nt = ~dE • E - ds • B, (2) 

where d£(a) is the electric (magnetic) dipole moment 
operator. Assuming the fields arc constant in time 
over the neighborhood of the atom, the cavity mode 
can be quantized [1] with E cx (a -|- a^)sin(^z), B cx 
i{a^ — a) cos(^2:), and H field = wa^a (with a and a) de- 
noting creation and annihilation operators of the mode) . 
Switching sin -(->■ cos by introducing b = ae~^^ , discretiz- 
ing the z-axis over the TV position states of the atom 
(^z|n) = '^j^\n)), and relegating the coupling strengths 
to a parameter A obtains the A'^-state Rabi Hamiltonian 




N-l 

H = uh'^h + A E (&e'T^" 6^e-*T^")|n)(n| 

rt=0 
N-l 

+ E E •^'^ + + I" + (3) 

k n=0 

For the two-state case (A^ = 2), this simplifies to the 
original Rabi model 

H2^Lob'^b + \{b + b^)a^ + Ja^, (4) 

where J = J n_ and ^ are the usual Pauli matrices. 
One can also interpret as a normal mode smeared 
over a tunneling A^-site system (discussed later). As a 
result, this extension maintains the correspondence be- 
tween atom-field interaction and a more general spin- 
boson model. 

It will now be shown that re-expressing the model in 
terms of the generalized spin matrices, the unitary gen- 
eralization of the Pauli matrices [26| . will reduce mathe- 
matical complexity while increasing physical understand- 



ing. Suppressing dependence on A^, generalized spin ma- 
trices for < j, k < N are defined (via modulo A^) as 

Sj,k - E + k\^ {Si,oy {So,it . (5) 

n=0 

With the details relegated to [2^, the reader need only 
keep in mind the function of the two indices: j determines 
the phase at each entry n while k determines the entry's 
location. The matrices 5i^o a-^d ^1 q elegantly express 
Hint while So^k + Sq ^ describes the neighbor couplings 
of Hsys- For < fc < K = [i (Af - 1)J (with [A^J the 
floor function), Eq. ([3]) is thus re-expressed as 

H = Ljbh + X{b S^^o + b^ Slo) 

+J5o,f +E'^fc(V + 4fe)- (6) 

k=l 

Energy representation. — One can now transform H 
into the energy representation H = HV using the uni- 
tary transformation 

N-l 

^ = 7^Ee'^''^'=.'=^i'0' (7) 

^ k=0 

linking the above formulation with the well-established 
picture of dipole transitions in A-state systems Q. The 
transformed Hamiltonian 

H = ujbh + X{bS2,ie"^ + b^f Sl^e-'^) + Hsys (8) 

models an A^-state atom coupled to a held mode (with the 
option for more modes jl^). The symmetry of the cou- 
pling determines which states are coupled by the mode 
and the state energies are determined by Hsys- For 
N = 2, Eq. ^ reduces to H2; the Rabi Hamiltonian 
in the energy representation [ax -H- CTz in Eq. ([4])]. The 
three- and four-state cases are reviewed below. 

A^ = 3: Setting Ji = AT in Eq. the three-state case 
in the energy representation is 

/iub^b + 2K Xbe-'^ A6t \ 
H3 = A6l-e^¥ Lob^b - K Xbe'^ . (9) 
V Xb Xbh^'^ Lub^-Kj 

The above is a three-level atom with an initially degen- 
erate ground state and energy separation 3K coupled 
to one cavity mode [2^. The states are thus arranged 
in a A conflguration (with inversion of K obtaining a 
V configuration), similar to well-studied A-systems [la |. 
However, dipole transitions occur between all three lev- 
els while extending the symmetry and maintaining the 
relative simplicity of the original Rabi model. The bot- 
tom left entry in Eq. ^ describes the process in which 
the atom makes a transition from the upper to the lower 



2 



level and a photon is annihilated Q • The RWA (with re- 
spect to ujb^b+Hsys) removes this transition, relating Eq. 
(|9]) to well-established extensions of the J-C model [l6[. 
The coupling between the ground states represents an 
ac-Stark shift (similar to the J-C model in the dispersive 
regime [soj . relevant to non-demolition measurements), 
which interestingly remains relevant after the RWA. 

N = A: For the four-state case, Hgys — 
diag {J -I- 2K, —J,J— 2K, — J}. Depending on the re- 
lation between J > and K [Fig. ffljc)], one can ob- 
tain either a double-A, tripod, or four-state configura- 

mill 



tion [191 131lj . Inversion of J obtains inverted tripod and 
double-A configurations; inversion of K leaves the system 
invariant just like inversion of J for N = 2. The cavity 
frequency w can be tuned to the three possible transi- 
tion frequencies of the atom, producing a four-parameter 
model for treating single- and (in the case) multi-level 
transitions in several related systems. Additionally, 
can be separated into two effective two-state systems as 
K — >■ 0. One striking feature is that the ground state 
can change for increasing values of A, a property not 
seen at < 4. Shown in Fig. EJb) for a particular 0- 
configuration, the ground state at small coupling (blue) is 
surpassed by the unperturbed first excited state (green) 
as the coupling increases. 

Conserved quantum numbers. — The generalized spin 
matrices allow one to easily construct the complete set 
of conserved quantum numbers for H, providing impor- 
tant insight into dynamics (§] and integrability Q. It 
can be shown that the Hamiltonian ([6]) possesses an A^- 
fold rotational symmetry and commutes with the rota- 
tions {T^nSo^n}nZo ' where the bosonic rotation TZn = 
exp {i'^nb^b) and the parity /reflection TZ = TZn is 
present for even A'^. These can be compiled into the gen- 
eral A^-state commuting operator 

K 

N - jTzs.^N + ^ Mn,So,k + nisi,), (lo) 
fe=i 

consisting of the family of k commuting Hermitian oper- 
ators multiplied by site couplings Jj. (with the additional 
parity operator for even A^). For the original N = 2 
case, N reduces to J multiplied by the well-known spin- 
boson parity axT^ 0, 0- This result shows that these 
A-state atom-cavity systems not only preserve parity for 
any even A^, but are classified by other quantum numbers 
for N > 2. For example, the three-state case contains a 
conserved quantum number S = 2,-1 while the four- 
state state system has two: parity p — ±1 and "cascade" 
number c = 0, ±2. 

Analytical insight. — The rotational symmetry of H 
allows decomposition into A^ infinite-dimensional sub- 
spaces {boson chains, denoted as i?Ar.„) via a group- 



theoretic transformation U l28[. In this partially diag- 
onalized representation, the Hamiltonian H = WHU is 



(a) 




0.2 



0.4 



0.6 



(b) 



0.8 

X/a> 



1.0 



1.2 



1.4 



diagonal in the spin subspace with {n'\H\n) 



Na 



FIG. 2. (color online) Correlation diagrams of energy 
vs. coupling A for (a) Hi, a resonant A-configuration with 
K = a;/3, and (b) H4, a 0-configuration with K — cu/A and 
J — ij/G {uj = 1). The numerical energies belong to chains 
HN,n (for n < N = 3,4) represented as red, blue, green, 
and cyan, in that order, while approximate S-GRWA ener- 
gies are dashed. The model predicts that the 0-configuration 
may have a dilTerent ground state at A/oj ~ 1 than at weak 
coupling. 



These chains are isomorphic to H and provide signif- 
icant numerical advantages [7|. For the two-site case, 
N — > Jctz, resulting in parity chains shown in the 
left panel of Fig. [TJa). The chains and their respective 
quantum numbers for the three- and four-state cases are 
depicted in Fig. W[h) and (c), respectively. The numeri- 
cal energies for ^3 and H4 are plotted in Fig. OJa) and 
(b), respectively, and each chain is labeled by a color. 
The spectrum of H demonstrates the familiar braid-like 
crossing pattern of the two-level Rabi Hamiltonian with 
the addition of more braids. 

The group theoretical approach is also useful for ex- 
tending analytical approximations, such as the symmet- 
ric generalized RWA (S-GRWA 0, applied in [ll). The 
S-GRWA energies (dashed in Fig. ^ are most accurate 
in the deep-strong coupling regime (A ~ uj), where the 
symmetry and chain structure are important While 
the S-GRWA also fares well at A ^ w, the symmetry is 
not terribly relevant in that region and the original RWA 
may be applied without loss of accuracy. In the weak cou- 
pling regime, it is anticipated that classical J-C collapses 
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and revivals [8| will occur in the system dynamics, but 
this time between multiple atomic states. The notably 
different behavior of H in the strong-coupling regime will 
likely be an extension of that described in [9| and will be 
chain-dependent. Both of these regimes (as well as tran- 
sitions between them) reveal opportunities for interesting 
manifestations of both well-known and newly-discovered 
phenomena of the original two-level case. 
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FIG. 3. (color online) In a manner analogous to [2l|], (a) is a 
pictorial representation of a symmetric three- well potential in 
the three-state limit, (b) is a visualization of the B800 ring in 
LH-II [s^l, representing a specific application of the A'^-state 
spin-boson to modeling multi-chromophoric energy transfer 
in periodic systems. 



site systems as previous approaches [2l|, [36| have revealed 
in the simplest two-site case. 

As a final note, instead of extending the number of 
modes (or even reservoirs [s^l), the A^-state Rabi model 
can be extended to many A'-state systems. This ap- 
proach would be similar to previous extensions of the 
two-state case [s^, but would include odd A^, potentially 
revealing phase transitions and other interesting physics. 

Summary. — This work introduces an extension of the 
two-state Rabi model [l| to describe dynamics of a more 
general A'-state periodic system. The symmetry of the 
system is utilized in a group-theoretical approach, re- 
vealing insight into its energies and conserved quanti- 
ties while also simplifying numerical analysis. A re- 
cently developed class of matrices [2^ provides an ele- 
gant method for obtaining the above results. Finally, the 
proposed infinite-mode extension generalizes the two-site 
spin-boson model 
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to dissipativc periodic A^-site sys- 
tems. 

Discussions with M. H. Devoret, S. M. Girvin, F. 
lachello, G. D. Scholes, and A. Nazir are gratefully ac- 
knowledged. I thank J. I. Vayrynen, K. A. Velizhanin, 
and D. Bokhan for help with preparation of this 
manuscript. This work is supported by an NSF Grad- 
uate Research Fellowship. 



Extensions. — Having examined the single- mode case, 
the dissipative version is now defined. As an extension 
of Leggett et al. 
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consider a continuous A^-well sys- 
tem with symmetric potential V{X^P) where the dy- 
namics is restricted to the A^-dimensional subspace of 
the well ground states. One then obtains Hgys by intro- 
ducing tunneling matrix elements Jk between the wells 
[see Fig. [3lja) for A^ = 3]. With the dissipative environ- 
ment approximated by a continuum of modes {qiTpi}, the 
A^-state spin-boson Hamiltonian is simply Eq. ([6]) with 
{&, w. A} — > {hi,uji, A;}. The interaction term satisfies the 
criteria of [33[ and simplifies to the degenerate two-site 
spin-boson model at A^ = 2. Other continuum normal 
modes can be added in a similar fashion j29| . 

Since the A^-state model preserves rotational symme- 
try, the dissipative H is an effective model for the sin- 
gle excitation manifold of a molecular ring or periodic 
chain interacting with a normal mode of a collective un- 
corrclated vibrational bath [23|. This model specifically 



includes the geometrical structure of the system, an im- 
portant property in excitation energy transfer [3J| . Cou- 
plings Jfc between all sites in the ring are included, al- 
lowing one to model systems with interactions other than 
nearest-neighbor. This version can model photoexcita- 
tion dynamics of molecular trimers (23 . 35 1 and larger 
rings. A specific exam ple is the 8-9 member B800 ring of 
photosynthetic LH-II [32|, illustrated in Fig. [2Ib). Re- 
cently developed methods [l5| for spin-boson dynamics 
can readily be applied to reveal similar insight into many- 
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Generalized spin matrices 

Suppressing dependence on N ^ generalized spin matrices [l| for < j, k < N are once again defined (via modulo 
N) as 

JV-l 

= ^ e'^^"^|n)(n + k\ = {Si^oY {So,if ■ 

n=0 

The matrices 6*1^0 said Sq^i are the generators of this unitary set and 5o,o is the identity. These matrices are ideal for 
describing periodic A^-level systems where the site interactions are circulant. While Sjq^ ■S'-j.o and 1^ = ^o.-fc, 

for arbitrary matrices 5*1^. = e^^^^ S-j_-k- For N ~ 2, Sq^i ~ (j^, S'1.0 = <^z, and Si^i = idy. For iV = 3, the relevant 
matrices are 



5o,i = 1 , So 
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Si,o = 





e 3 
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52,0 = 





- 271 












e 3 




S,i = e'^ , 5i,2 = 




^2,1 =00 e-'- , 52,2 = 



Letting a, b,j, k be integers, one obtains a phase when two matrices are permuted: Sa.bSj^k = e* « ''"''^ Sj^kSa,b- An 
additional identity useful in manipulating these matrices and proving the steps below is 



N-l 

1 ^ e^¥"('-™) = Si,^. 



n=0 



Partial diagonalization of the A'^-state Rabi model 

Since the A^-state Rabi model contains an iV-fold rotational symmetry, one can partially diagonalizc H in the site 
subspace using the prescription of 0]. The generalized spin matrices significantly simplify this prescription and the 
general transformation is a composition of spin and boson operators: 

JV-l 

^ to 

This unitary U commutes with the coupling term Hint since the phase contribution from permuting the spin matrices 
exactly cancels the phase from permuting the boson operators. The transformation partially diagonalizes the system 
matrix Hgys since U^So^kU = TJ.j^S'^g. The transformed H is diagonal in the spin space with {n'\H\n) = Sn',nHN^n 
and boson chains 

K 

HN,n = {n\U^HU\n) = ujb'^b + A(6e*T?" + 6^6"'^^") + (-1)" J7^ + Jki'R-ke'^"'' + 7^^e-*^?"'=). 

k=l 

The transformed form of the commuting operator N is simply the diagonal form of Hgyg , namely, 

N = Hsys = JSn ^ + MSk,o + Slo)- 



For the nth diagonal of N, the^ coefficients of the site couplings Jk arc conserved quantum numbers that de- 
scribe the corresponding chain H^^n- As an example, N = Juz when N = 2, revealing chains of parity ±1. 
Likewise, N = diag{2/v, — A', — A'} for iV = 3, revealing a conserved quantum number 5 = —2,1. For = 4, 
N — diag{ J + 2 A', — J, J — 2K, — J} and there arc two c^uantum numbers: the parity p — ± and the cascade number 
5 = 0, ±2. The former determines whether the state is in the degenerate — J chain while the latter resolves the 
remaining two chains J±2A'. Note that there are no quantum numbers which can be used to distinguish the negative 
parity chains for the four-state case. 



Symmetric generalized rotating wave approximation (S-GRWA) 

The 1-by-l S-GRWA [i| is derived by applying the generalized displacement operator 

2?(5i,oA/c^) = cxp[(65i,o - b^Sl„)X/uj] 

to H and writing out each transformed chain of HT> in the boson Fock space. This removes the linear bosonic 
coupling terms and (keeping only diagonal elements) obtains energies (for m = 0, 1, 2, . . .) 



= - ^ + Je-2^^/-^L„,(4AV^')(-l)"+'" + jZ 2 J^e-^l-^/^l' A,„ cos (^fc (n + m) + 6,) 

^ k=l ^ ^ 

where the third term is only for even A^, Lm{x) is a Lagucrre polynomial, and 

Ok = {X/ujfsm{27rk/N) 
ak ^ I - e ' . 

For N = 2, the 1-by-l S-GRWA simplifies to well-known results (see [3l and refs. therein). Interpreting as a 
tunneling A"-site system in the position representation, the energies consist of contributions from standing waves 
along the ring/chain of sites weighted by site coupling Jk, displaced by 9k, and decaying as O(A^). The corresponding 
eigenfunctions consist of the position basis rotated along the ring by W and displaced vibrationally by I?^ . While the 
opposite-site (parity) interaction J is the only interaction in the two-state case, it naturally becomes less relevant as 
2N becomes large. The two limits of the system are A ^ {Jk}, which simplifies H into a collection of symmetrically 
displaced harmonic oscillators, and A <^ {Jk}, which results in an uncoupled A^-state system and an oscillator. 
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